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The known correspondence between the Kronig-Penney model and cer- 
tain Jacobi matrices is extended to a wide class of Schrodinger opera- 
tors on graphs. Examples include rectangular lattices with and without 
a magnetic field, or comb-shaped graphs leading to a Maryland-type 
model. 



I Introduction 

Schrodinger operators on L'^iT) , where F is a graph, were introduced into quantum 
mechanics long time ago [|l[] . In recent years we have wittnessed of a renewed interest 
to them — see [0, |^, ^, |, |^, |^, ^ |^, |10|, |11|, |12[ and other, often nonrigorous, studies 
quoted in these papers — motivated mostly by the fact that they provide a natural, 
if idealized, model of semiconductor "quantum wire" structures. 

Jacobi matrices, on the other hand, attracted a lot of attention in the last 
decade, in particular, as a laboratory for random and almost periodic systems. The 
most popular examples are the Harper and related almost Mathieu equation dating 
back to HTsj [14|, [T5|] ; for more recent results and an extensive bibliography see, e.g., 
T6|, O, [181 O, |20, pl|] . The underlying lattices are mostly periodic of dimension one 



or two; however, more complicated examples have also been studied p2|, 

In case when F is a line with an array of point interactions, i.e., the Schrodinger 
operator in question is a Kronig-Penney-type Hamiltonian, there is a bijective cor- 
respondence — dubbed the French connection by B. Simon [TB 
systems and certain Jacobi matrices 



between such 



18, 24, m 



28]. The aim of this paper 



is to show that the same duality can be established for a wide class of Schrodinger 
operators on graphs, including the case of a nonempty boundary. In general, the 
resulting Jacobi matrices exhibit a varying "mass" . 
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II Preliminaries 



Let r be a connected graph consisting of at most countable families of vertices 
V = {Xj : j e 1} and links (edges) C = {Cjn ■ {j,n) G Ic C / x /} . We 
suppose that each pair of vertices is connected by not more than one link. The set 
Af{Xj) = {Xn : n G z/(j) C / \ {j}} of neighbors of Xj , i.e., the vertices connected 
with Xj by a single link, is nonempty by hypothesis. The graph boundary B consists 
of vertices having a single neighbor; it may be empty. We denote by and Jj 
the index subsets in I corresponding to B and the graph interior X := V \ B , 
respectively. 

r has a natural ordering by inclusions between vertex subsets. We assume that 
it has a local metric structure in the sense that each link Cjn is isometric with a line 
segment [0, £jn] ■ The graph can be also equipped with a global metric, for instance, 
if it is identified with a subset of M'^. Of course, the two metrics may differ at 
a single link. Using the local metric, we are able to introduce the Hilbert space 
L^(r) := 0(j>)e/£ -^^(0, ijn) ] its elements will be written as ip = {ipjn '■ (j, n) G Ic} 
or simply as {ipjn} ■ Given a family of functions V := {Vjn} with Vjn G L°°(0, ijn) , 
we define the operator He = H(T,C, V) by 

Hci^Pjn} := { -^p-n + Vjn^.n : (j, u) E Ic} (2.1) 

with the domain consisting of all ip with ipjn G W'^''^{0,£jn) subject to a set C of 
boundary conditions at the vertices which connect the boundary values 

V'in(j) := lim^tpjuix) , Vin(j) := }}^_^'^jnix) ; (2-2) 

we identify here the point x = with Xj . 

There are many ways how to make the operator self-adjoint; by standard 



results [0, ^ each vertex Xj may support an Nj parameter family of boundary 
conditions, where Nj := cardz/(j) . These self-adjoint extensions were discussed in 
detail in |^; here we restrict ourselves to one of the following two possibilities which 
represent in a sense extreme cases [jlOl , |ll]: 



(a) 6 coupling: at any Xj G X we have 4'jn{j) = ^jm(j) =: "^j for all n,m & ^{j) , 
and 

Y: ^',nij) = (2.3) 
n&u{j) 



for some aj E IR. 

(b) (5; coupling: Vjn(j) = V'jm(j) =: ^'j for all n,m G z/(j) , and 

Y: ^,n(j) = (2.4) 

n&u{j) 

for pj eR. 

The relations ( |2.3| ) and ( ^.4] ) are independent of V , since potentials are supposed 
to be essentially bounded. At the graph boundary we employ the usual conditions, 

ipj cos LJj + ip'j sin Uj = 0, (2.5) 
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which can be written in either form (infinite values allowed); for the sake of brevity 
we denote the set C as a = {aj} or P = {Pj} , respectively, with the index running 
over V, and use the shorthands Ha ■= H(T, {aj},V) and Hp := H(T, {Pj},V) . 
Looking for solutions to the equations 

Hctp = C = a,P, (2.6) 

we shall consider the class Dioc{Hc) which is the subset in \/{j,n)£ic -^^(0,^jn) (the 
direct sum) consisting of the functions which satisfy all requirements imposed at 
ip G D{Hc) except the global square integrability. 

The conditions ( [^.3D and ( [^.4D define self-adjoint operators also if the coupling 
constants are formally put equal to infinity. We exclude this possibility, which 
corresponds, respectively, to the Dirichlet and Neumann decoupling of the operator 
at Xj turning the vertex effectively into Nj points of the boundary. On the other 
hand, we need it to state the result. Let H^ and H^ be the operators obtained 
from Ha, Hp by changing the conditions (|2.3|), ( p.4| ) at the points of X to Dirichlet 
and Neumann, respectively, while at the boundary they are kept fixed. We define 
fCa '■= {k : E a{H^), Imfc > } and ICp in a similar way. 

Ill Main result 

Consider the operators Ha, Hp defined above. We shall adopt the following as- 
sumptions: 

(i) There is C > such that ||Vj„||oo < C for all (j, n) G Ic- 

(ii) io := inf { ijn : (j, n) e Ic} > . 
(lit) Lq := sup{ ijn ■ (j, n) e Ic} < oo . 
(iv) Nq := max{ cardz/(j) : j G / } < oo . 

On Cnj = [0, ijn] (the right endpoint identified with Xj ) we shall denote as 
Mj„, C = a,P, the solutions to — /" + Vjnf = k"^ f which satisfy the boundary 
conditions 

U%{^,n) = l-«J'(^,n) = , !-<(£,.) = {u%)'{^,n) = 

and 

t;;„(0) = l-(t;p'(0)=0, l-t;f„(0) = (<)'(0)=0 if n G Jx , 
w^„(0) = sinu;„ , {v%y{^) = - cosu;„ if neh] 

their Wronskians are 

Wf^ = -<„(£,„) , = (<)'(£,„) , 
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respectively, or 

W^;n = <n(0), l< = -(<)'(0) if nG/x, 

Wf„ = -mJ„(0) coscun - (mJ„)'(0) sinci;^ if n G • 

If not necessary we do not mark explicitly the dependence of these quantities on k . 

Theorem III.l (a) Let ip G Dioc{Ha) solve ( \2.(^ for some k ^ fC^ with k"^ G IR, 
Im A; > . Then the corresponding boundary values satisfy the equation 

neu{j)nlx \n£u{j) jn J 

Conversely, any solution {ipj : j G Ij] to ( ^Tlj ) determines a solution of ( j^T^j by 
^U^) = ^ u%{x) - ^ vf^ix) if n G n h , (3.8) 

jn jn 



^.n(^) = -^<n(^) if nEuU)nls. (3.9) 



(b) For a solution ifj G DiodHp) of ( ^.(p with k ^ JCjj , the above formulae are 
replaced by 

E S + f E ^ + a]*; = o. (3.10) 



and 



^,n(x) = -^u%{x) + ^vlix) if n G n Jx, (3.11) 

''^ jn jn 



i^jnix) = ^ v^^ix) if n G z/(j) n Is . (3.12) 



(^cj Under (i), (ii), ip ^ L'^iX) implies that the solution {'ipj}, {i^j} of p. 1[ ) and 

■!(] ), respectively, belongs to ■ 
(d) The opposite implication is valid provided (Hi), (iv) also hold, and k has a 
positive distance from from Kq ■ 

Proof. (a,b) We shall consider throughout, the argument for Hp is analogous; 
for simplicity we drop the superscript a . If n E Ij , the transfer matrix on Cnj is 



T„,(x,0) = W-. 



-1 

jn 



( Ujn{x) - u'j^{0)Vjn{x) Ujn{0)Vjn{x) 
V U'j^ix) - U'j^{0)v'j^{x) Ujn{0)v'j^{x) 
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the Wronskian is nonzero for k ^ /Cq . This yields an expression of ipjn{x) in terms 
of ipjn{n) =: ipn and ipjnin) , in particular, 

... l-^jn(0)^jn(^i«) , ^ / WW ^ 

-'iPjnU) = + V(^i'^)^ir^N ; 

the sign change at the Ihs of the last condition reflects the fact that (|2.2| ) defines the 
outward derivative at . We express ipjn{n) from the first relation and substitute 
to the second one; this yields 



W~n 

If n E Ib we have instead 



3 ' 



VjnU) - TJT Vn + ... WjnK''^) > 

we may write, of course, ipl^ := ipjn{n) . Using the relations ipn cos Un+ip'n sinc<j„ = 
we find ^ 

so ( p.7| ) follows from ( p.2|) . The transfer-matrix expression of ipjnix) together with 
the mentioned formula for i/jj^in) yield (p.8|); ( ^l9|) is checked in the same way. 
(c) Without loss of generality we may consider real ip only. Any solution to the 
Schrodinger equation on Cnj can also be written as 

where Wjn is the normalized Neumann solution at Xn , Wj„(0) = 1 — Wjn{0) = 0. 
The argument of |]T^ Sec. III. 2.1] cannot be used here, even in the classically allowed 
region. Instead we employ a standard result of the Sturm-Liouville theory pO| , 
Sec. 1.2] by which 

sinkx sin ky 

Vjn{x) = —j^ + KD{x,y)—^dy, (3.13) 

and similar representations are valid for wjn and their derivatives; the kernels have 
an explicit bound in terms of the potential Vjn ■ The latter is essentially bounded, 
by (i) uniformly over C . Hence there is a positive ii < ^io such that 



max 



{\Vjn{x)\, \v'j^{x)-l\, \Wjn{x)-l\, \w'j„{x)\} < ^ 



for X G [0,ii) and any {j,n) & Ic] we infer that 
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An analogous estimate can be made for ijjj^{n) ; summing both of them we get 

li ijj e L^{T) , it belongs to D{H^) , so ^Ij' G L'^iV) also holds. Let Cjn{.^i) be the 
link Cjn with the middle part {(!-i,(-jn — ^i) deleted; then 



2 

L2(r)- 



^'Ili2(r) 



E 



> 



> 



25 



4£i 



(i,n)6-f£ 
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which yields the result. 

(d) First we need to show that has a uniform bound for k ^ JCa Without loss 
of generality we may suppose that inf /Cq, > ; otherwise we shift all potentials by 



a constant. In view of the Sturm comparison theorem [31, Sec.V.6], see also [29 
Sees. XIII. 1,15], the zeros of Vjn{ijn,-) and Vj^{£jn,-) form a switching sequence; 
hence it is sufficient to prove that '^^^"^^^ is bounded uniformly away from zero 



around each root ko G /Cq of Wjnik) = 0. Using the known identity []T9|, Sec. 1. 5], 

dv'j^{x,k) 



29|, Sec.XIII.3], 



Vj„{x, k) 



at X 



dk^ ^ 
^-nj, k = ko , we find 



dk^ 



1) ^ doc . 



dVjnjijn, kp) 

dk 



2fco 



'^jni.^jnj koj J Cjn 



Vjn{x, koY dx . 



We have supposed that ko > inf /Cq, > ; furthermore, by the argument of the 
preceding part we have J^^^Vjn{x, ko^dx > -^il- At the same time, using the 
representation ( |3.13| ) we find that \Vj^{ijn, ko)\ < Ci for a positive Ci independent 
of j, n . Hence to a given k ^ /Cq there is a C2 > , again independent of j, n , 
such that Wjn > C2 ■ The representation ( |3.13|) also yields another uniform bound, 

C3 := sup I [[nj„[[L2, [[i;jn||L2, : (j, ^ Ic] < 00. 

6 



The relations (|3.8|) and ( p.9|) together with (iv) now yield 

" = E U.nV < E \W,n\-'{m'\\u,r^?+\n^\\v,nf} 



SO the result follows. | 

Remarks. 1. Since the operators He are below bounded, one can express for 
k"^ e p{Hc) n p{Hq) the resolvent difference {Hq- k'^)~^- {Hq- where Hq = 

H^, Hg , by Krein's formula [p^, |3^, this expression becomes singular under 



'/3 

the condition ( |3.7|) or (|3.1CI|) , respectively. 

2. If two vertices are joined by more than a single link, the theorem can be applied 
if a vertex with the free 6 coupling, aj = , is added at each extra link. Vertices 
with just two neighbors are also useful, if we want to amend the potentials Vjn with 
a family of point interactions — for an example see ^5 



3. If (i)-(iv) are relaxed, the imphcations (b), (c) may still hold with replaced 
by a weighted C"^ space. If V is equipped with a global metric, one can establish a 
relation between the exponential decay of ip and of the boundary- value sequences. 



IV Examples 

Example IV. 1 (rectangular lattices): Let F be a planar lattice graph whose basic 
cell is a rectangle of sides j = 1,2, and suppose that apart of the junctions, the 
motion on F is free, Vjn = . By the proved theorem, the equation (2^) leads in 
this case to 



■'n,m+l n+l,m + 4>n-i,m) sin ki2 



(4.14) 



(l>nm = 

in the 5 and 5'^ case, respectively, where 

:= - ^ sin Hi sin H2 , V^^{k) := - P^mk sin Hi sin H2 , (4.15) 
k 

and (j)nm = i^nm, i^'nm ■ particular, for a square lattice we get 

{hQ^)nm - COS k£ + sin k£^ tpnm = 

and an analogous equation in the 6^ case, where ho is the conventional two- 
dimensional Laplacian on . On the other hand, if ii ^ £2 , the free operator 
has a periodically modulated "mass" which leads to nontrivial spectral properties 
even if the underlying graph Hamiltonian is periodic, i.e., anm = a or Pnm = P , 
with a non-zero coupling parameter. A detailed discussion of this situation can be 
found in [0, |ll|, |12|. 
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Example IV. 2 (magnetic field added): Suppose that F is embedded into M'^ in 
which there is a magnetic field described by a vector potential A . The boundary 
conditions (|2.3|) and ( p.4|) are modified replacing ipjnij) by ipjnij) +iAjn{j) , where 
Ajn{j) is the tangent component of A to £j„ at Xj ; we suppose conventionally 
that e = —1 . For the 6 coupling this is well known 0; in the 6^ case one can 
check it easily. 

There in no need to repeat the above argument, however, since the magnetic 
case can be handled with the help of the unitary operator U : i^^(F) -^^(r) 
defined by 



X 



where Ajn is again the tangent component of the vector potential and fixed 
reference points. Then the functions {Uip)jn satisfy ( p.3| ) and (|2.4|), respectively, 
and the equations ( |3.7[ ) and ( ^.1U| ) are replaced by 

provided the magnetic phase factors Aj obey the consistency conditions 

Aj- An = / Ajn{y) dy 



following from the continuity requirement. 

In particular, if F is a rectangular lattice considered above, Anm '■= ^{n—m)/2 
corresponds to a homogeneous magnetic field in the circular gauge, with the fiux 
$ := Biii2 through a cell, and no other potential is present, the equation ( [4.14 ) is 
replaced by 

(4.16) 

+ {vLik) T 2 sin k{ii+ £2)) 



where the discrete potential is again given by ([4.15|) . If = £2 and the coupling 



constants vanish, anm = f^nm = , we obtain in this way the discrete magnetic 
Laplacian of (or the Harper operator of [|13|, |19| for the Landau gauge), while for 



a general rectangle the parameter A representing the anisotropy in |^ is replaced 
again by the periodically modulated "mass" . 
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Example IV. 3 (comb-shaped graphs): Let F consist of a line to which at the 
points Xj := jL line segments of lengths ij are joined; at their ends we impose the 
conditions (|2.5| ). The equation ( p.6|) now yields 

{ho<j))j + (Vf{k) T 2 cos kl) <f)j = (4.17) 

in the 6 and 6'^ case, respectively, where 0j = ipj, ipj , and 

In particular, in the absence of an external potential on F the last relations become 

= - (^cot{kij+r]j) + ^) sin kL , Vf{k) = - {ta.n{kij+7]j) + Pjk) sin kL , 

where rjj := arctan(/c tantUj) . If the coupling is ideal, aj = or Pj = , the loose 
ends correspond to the Dirichlet condition, uj = , and the "tooth" lengths are 
ij := \j\i , we get thus an equation reminiscent of the Maryland model ^] with 
a fixed coupling strength and an additional periodic modulation of the potential. 
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